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Abstract 


The  position  of  a point  on  a surface,  as  determined  by 
measurement,  is  subject  to  errors  of  observation.  The  direction 
and  distance  of  the  measured  point  from  the  true  point  depends 
not  only  on  the  errors  of  the  measured  coordinates  but  also  on 
the  properties  of  the  two  families  of  lines  of  which  the 
measured  coordinates  are  the  independent  parameters.  Practical 
considerations  generally  preclude  the  use  of  measured  coordinates 
for  which  the  error  of  position  depends  only  on  the  errors  of 
observation,  or  for  which  the  effect  of  the  coordinate  system 
is  easily  understood. 

When  the  distribution  functions  of  the  errors  are  known, 
it  is  possible  to  construct  around  any  point  representing  the 
true  point  on  tne  surface  a closed  curve  along  which  the 
probability  of  occurrence  of  a corresponding  measured  point  is 
uniform,  and  inside  which  there  is  a specified  probability 
that  a measured  point  will  fall.  Accordingly,  the  same 
closed  curve  may  be  drawn  around  a measured  point  to 
represent  the  aize  and  shape  of  the  smallest  region  inside 
which,  with  the  specified  probability,  the  (unknown)  true 
point  will  lie. 

Details  are  worked  out  for  two  kinds  of  error  functions  -- 
the  normal  law,  and  errors  of  limited  size.  For  the  former. 


the  paper  shows  how  to  calculate,  for  each  position  of  a 
measured  point  on  the  surface,  the  boundary  inside  which  there 
is  a specified  probability  that  the  true  point  will  lie.  In 
the  latter  case,  equations  are  given  which  define  the  smallest 
region  inside  which  the  true  point  will  certainly  be  found. 

While  the  analysis  has  been  restricted  to  the  plane, 
the  ideas  are  applicable  to  other  surfaces  and  to  three 
dimensions . 
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Inti  -c.  action 

'he  coordinate  system  employed  to  find  the  position  of 
a point  by  physical  measurement  is  usually  determined  by  necessity, 
convenience,  or  accuracy  in  making  the  measurement.  The  choice 
must  often  be  mathematically  inconvenient.  In  two  dimensions 
each  measured  coordinate  is  a parameter  of  one  of  two  families 
of  loci.  A pair  of  coordinates  determine  a particular  locus 
from  each  of  the  two  different  families  and  these  two  curves 
intei-t.ee t to  determine  a point. 

Consider  the  two  coordinate  curves  which  pass  through  the 
actual  point  and  another  pair  of  curves  whose  parameters  are 
erroneous  measured  coordinates . These  four  curves  form  what 
may  be  called  a curvilinear  parallelogram,  with  the  point  at 
one  vertex  and  the  erroneous  point  at  the  opposite  vertex. 

The  distance  between  these  two  points  dependn  on  the  sides  of 
the  parallelogram  and  the  angles  ?t  its  vertices.  These 
elements  in  turn  depend  on  the  size  and  relative  signs  of  the 
errors,  the  nature  of  the  coordinate  system,  and  the  values  of 
the  coordinates  at  the  point.  The  re3atlons  between  an  error 
of  petition  and  the  errors  of  observation  from  which  it  arises 
are  much  more  complicated  In  general  than  they  are  in  the 
cart:; oian  system  or  even  in  orthogonal  curvilinear  systems. 

Since  it  Is  often  Impossible  to  avoid  using  complicated 
coordinate  systems , it  is  necessary  to  understand  these  matters 
if  ~ ? is  to  make  u choice  of  measuring  tools  and  their 
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locations  which  will  satisfy  the  requirements. 


Probability  considerations 


Let  U and  V be  the  curvilinear  coordinates  of  a point 
in  the  X,Y  plane;  and  let  TJ  + v,  V + v be  erroneous  measured 
values  of  the  same  pair  of  coordinates.  Assume  that  U and  V 
have  the  same  physical  dimensions.  Let  their  errors,  u and  v, 
be  of  normal  distribution  and  of  known  precision.  The 
probability  that  an  error  in  U lies  between  u and  u + du  while 
the  corresponding  error  in  V lies  between  v and  v + dv  1b 
then  the  function  , 


e 


/ L 


(1) 


where  the  k‘s  are  experimentally  determined  moduli  of 
precision.  The  set  of  all  points  for  which  this  probability 
is  uniform  are  subject  to  the  condition 


(2) 


The  parameter,  h,  may  be  determined  by  substitution  in  (2) 
of  any  particular  pair  of  values  of  u and  v.  The  corresponding 
value  of  the  probability,  (7.),  is 


ft*  — "iC 


(3) 
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and  this  value  is  uniform  over  the  band  between  the  two 
ellipses  wnoBe  parameters  are  h and  h + dh.  Since  the  area  of 

the  ellipse,  (2),  is 

-JT  K 

A = ~C 


the  area  of  the  region  between  the  ellipses  Is 


<L  A = -~r7  ^ 

1< 

Changing  the  element  of  area  in  (3)  from  dudv  to  dh.  we  find 
the  probability  of  a point  being  inside  any  particular  ellipse 
is 

p ~ = / ~ g 

c 

(4) 

If  k = Vc*  , the  ellipses,  (2),  become  circles 


y-  f •*-_  A,  \ 

iT  - ? C ^ =TZV  (5) 

for  which 

p=  I-  e 

r ( 6) 
Imagine  that  v?e  have  made  a very  large  number,  N,  of 
pairs  of  measurements  'U  + u , V + v)  of  the  coordinates  (U,V) 
of  a single  point.  Suppose  all  these  measured  points  to  be 
plotted  in  the  U,V  plane  and  also  in  the  X,Y  plane.  In  teimc 
of  the  N plotted  point;]  in  the  U„V  plane,  the  above  equations 
have  the  following  interpretation.  The  density  of  points  per 
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unit  area  is 


/v 


i 


/c 


Over  the  narrow  band  included  between  two  ellipses  of  the 
family. 


*2—  r,  4 *3  , /) 

jjl  +iL'  v~  ■=.'*-* 


whose  parameters  are  h and  h + dh,  the  density  of  points  is 

uniform  and  the  number  of  points  included  in  the  band  is 

Ne~  dh.  Suppose  we  divide  the  space  around  (U,V)  into  a number 


of  such  adjacent  elliptical  bands 


fcy  constructing  several 


concentric  ellipses  of  the  family  (2).  The  probability  that  a 
single  measured  point  will  fall  outside  one  of  these  ellipses  is 
a simple  function  of  the  parameter  h,  as  we  have  seen  in  (4). 

The  picture  drawn  here  in  the  U,V  plan®  is  exactly  the  same  as 
we  should  have  in  the  X,Y  plane,  had  we  been  measuring  X and  y 
directly  Instead  of  indirectly  through  U and  V which  are  functions 
of  X and  Y.  The  question  of  practical  interest  is:  what  does 

this  plot  in  the  U,V  plane  look  like  when  it  is  mapped  into  the 
X, Y plane?  Under  any  point  transformation  X « X(U,V},  Y = Y(U,V), 
we  should  expect  the  ellipses  (or  circles,  vjhen  k = k* ) to  map 
into  some  kind  of  closed  curves  with  the  Images  of  the  same 
points  as  before  lying  between  them.  In  other  words,  after 
the  mapping,  fne  probability  of  a single  measured  point  falling 
outside  a given  closed  curve  in  the  X,Y  plane  is  the  same  as 


that  of  its  failing  outside  the  image  of  this  c 


trm  ve  in  the  TJ , V 
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plane.  It  will  be  shown  that,  subject  to  very  mild 
restrictions,  small  ellipses  (or  circles)  in  one  plane  always 
reap  into  ellipses  in  the  other  plane.  The  size,  orientation, 
and  eccentricity  of  the  mapped  ellipse  generally  vary  from 
point  to  point  in  the  plane.  We  now  proceed  to  examine  the 
nature  of  these  transformations. 


Equations  for  the  mapping  of  small  regions 


Let  the  measured  coordinates  (U,V)  be  related  to  the 
cartesian  (X,Y)  through 

D-=.v(y,Y)  v = v/  Xj'f) 


or  the  inverse 


X* 


X(' 


, , \ 


V' - Y(uJ\,') 


(8) 


in  which  the  functions  are  all  supposed  to  be  real.  Either  (7) 
or  (8)  provide  all  necessary  information  for  mapping  one  plane 
into  the  ocher;  but  since  we  are  interested  only  in  mapping 
separate  small  regions  we  may  take  advantage  of  the  fact  that 
the  transformation  for  small  regions  is  generally  linear.  In 
place  of  (7)  or  (8),  which  may  be  of  any  form,  ve  may  substitute, 
linear  transformations,  the  elements  of  whose  .matrices  are  the 
partial  derivatives  of  (7)  or  (8)  evaluated  at  3orae  fixed  point 
in  the  small  region  being  mapped. 


If  X and  Y may  be  expanded  in  Taylor  series  in  the 
vicinity  of  a point  (U^  ) , v;e  have 
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Z-  - 


*t.e 
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X X _ i 

i FX_ 

- IL. 

bU 

+ Fv-  * 

( b'u2- 

>■  F u. 

TY  u 

*■  TO 

by  z\ 

\ 

b u*- 

LL~4-  X. 


i Li1-  4- 


> ? r„J~  nv  -h  - *. 


Fv** 

a’f  ,4) 


4- 


(9) 


In  which  x =»  X - Xq  , y « Y - YQ,  u .=  U - UQ,  v =■=  V - VQ  are 
finite  Increments  which  we  may  Identify  with  the  errors  under 
discussion.  In  (9)»  the  partial  derivatives  are  to  be 
evaluated  at  (UQ,V0)  . With  proper  restrictions  as  to 
continuity  and  size  of  the  region  we  may  therefore  employ  the 


equations 


t /■ 

j 


iFL  !l  + ir- 

> U c V 

X V.  u +.  ilX  u- 

hU  * V 


(10) 


and 


xc  — 


IT  — 


^ X 

"a  V 

> X 


* 4- 


IF 


X 


ay 
> v 
3Y^ 


(11) 


in  which  the  coefficients 


£ U 


, etc.  are  constants. 


In  order  to  shorten  the  notation  and  emphasize  the  linear 
character  of  the  approximations , we  re-write  (10)  and  (11)  as 

J / 

y - a-  lu  +.  &_  v — 

//  / 7L 

^ / i 

i =-  a,  .cc.  <r?  v— 

X/  2.Z. 


(12) 


_o_ 


and 


U.  = 

v--  ^ * + Ct^tr 


d» 


in  which  ail  the  <21...  are  real. 

V 

While  we  may  use  these  approximations  freely,  because 
they  are  generally  accurate,  it  is  true  that  one  could 
encounter  c ircumstancea  in  which  the  parameter  h is  so  large 
and  the  successive  partial  derivatives  of  such  values  at 
particular  points  of  the  field  that  (12)  and  (13)  would  be 
inaccurate.  However,  the  main  purport  of  what  we  have  to  say 
about  the  probabilities  associated  with  the  images  in  the  XPY 
plane  of ' the  ellipses  (2)  in  the  U„V  planes  13  unaffected  by 
whether  we  can  employ  (12)  and  (13)  or  must  use  (7)  and  (8). 

The  important  fact  is  that,  given  any  probability 
functions,  we  can  calculate  the  contours  of  constant  proba- 
bility in  the  U„V  plane  and  then  draw  the  corresponding  contours 
of  the  same  probabilities  in  the  X,Y  plane.  For  inonotonic 
decreasing  probability  functions,  the  latter  show  the  size 
and  shape  of  the  smallest  regions  in3ld.e  which  there  is  a 
specified  probability  that  a measured  point  w ill  occur.  It 
is  evident  that  this  can  always  be  done,  regardless  of  the 
shape  of  the  mono tonic  error  functions,  or  what  special 
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devlces  are  employed  to  effect  the  transformation.  Thus  the 
basic  ideas  of  this  discussion  are  independent  ox'  any 
particular  assumptions  such  as  (1)  or  the  use  of  linear 
transformations . 


Elements  of  the  ellipses  in  the  X,Y  plane 


The  ellipse  (2),  in  the  U,V  plane,  becomes  in  the  X,Y 


plane 


/ -/-/A  V* 

'fL  * L--  / ^ 


The  discriminant  of  the  quadratic  form  on  the  right  is 


which  shows  that  an  el-*  ipse  transforms  into  an  ellipse. 
Obviously  this  is  true  regardless  of  the  direction  of 
trans format ion . 

Let  us  write  (15)  ih  the  more  compact  form 


* Note  that  E,  F,  and  G as  here  defined  are  not  identical  with  the 
familiar  S,  P and  G of  differential  geometry  except  when 
k ® \C»  l. 
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and refer  it  to  new  coordinates  x* ,y'  so  that  it  takes  the 

form 


•7 


This  we  can  do  by  means  of  a rotation  about  the  point  Xq,y^  - 

A = k ' y l<u^>  S- 

y - S"  -+  y'‘  c-cru  ^ 


such  that 


1 

- s 


»9-  = - 


£.-&  ' 


( O '<  Z-&  '<  7H  ) 


As  E and  G are  positive,  we  find  in  the  relations 


<rvxj  = 


ztr  XL 


Z&  =r- 


ztXL 


t9" 


■( 


2 «. 


'} 


M 


K T (f-  trj)  V''- 

S 


ijrf'. 


^ /\J(£-if)  ' x-  ‘ILF  = -h  ys/(fv-6jZ-X-  O 


that  the  upper  sign  is  to  be  taken  when  F > 0 and  the  lower 
sign  when  ? < 0.  The  numbers  a and  b in  (ly)  are 
be  given  by 


U7) 


(18) 


U9) 


(20) 


found  to 
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where  the  choice  of  signs  depends  on  the  sign  of  F In  the 
same  way  a3  before. 


Vie  note  that 

R = 0 

implies  that 

E = G ; F « 0 

These  are  the  conditions  for  the  ellipse  to  be  a circle.  The 
eccentricity  of  the  ellipse  is 

'A 


o \ Y'l 

! i 

e = J 


If  the  conditions 


£ ¥=  & 


are  satisifed  the  locus  Is  an  ellipse  whose  axes  are  parallel 
to  the  axes  of  X and  Y and  whose  eccentricity  is 

e §")  > (*<*■>  j e = (/-  fr)  ^ (£  > 

The  major  and  minor  semi-axes  in  the  general  case  are 


cL 


£ +-  & + K 


v 'ir 


(21) 


(22) 


(23) 


or 


(e  < *) 


The  chief  results  of  this  section  are  the  expressions 
for  major  and  minor  serai-axes,  (23),  and  the  direction  of  the 
x*  axis,  (20).  The  rule  of  signs  implies  that  when 
F > 0 the  major  axis  lies  in  the  second  quadrant,  and  when 

F < 0 it  Is  in  the  first  quadrant. 


Bounded  errors 


Instead  of  following  the  normal  distribution,  it  may  be 
found  that  errors  of  absolute  value  greater  than  a certain 
small  number  do  not  occur.  For  example,  a set  of  measurements 
of  an  angle  by  an  observer  using  a certain  sextant,  may  show 
that  the  chance  of  an  error  greater  than  two  miru.-  ec  is 
negligible.  If  this  Is  true,  we  may  employ  as  the  smallest 
contour  of  unit  probability  in  the  U,V  plane  the  rectangle 
formed  by  the  four  lines 

U,  = U-o  Il.-z.-~ 


*tr  = VS 


'LT  - 


- -LT„ 


(24) 


which  maps  into  the  oblique  parallelogram  bounded  by  the  lines 


u,  - a.  *•**-,,.  y ) 


(25) 


The  coordinates  (x,y)  of  the  four  vertices  of  the 
parallelogram,  referred  to  (Xo,YQ)  as  origin,  are  given  b-y 
the  inverse  transformation,  (12), 


J > ( ic  - =±  ^ , 


Vc±Vo 


(26) 


and  the  squares  of  the  semi-diagonals  are  given  by  the  two 
values  assumed  by 
,2 


* 


s*  (a;>0  ^v- 

~ £.  U'  -t-Z-J’**-  V~  -t-'Zf  V-  , j u^Xr=  w.0V“  > 0 1 


(2?) 


( 


cLo-T-ci-  u.ir  = c.  o/ 


With  proper  attention  to  sign,  the  angles  at  the  vertices 
of  the  parallelogram  are  given  by 


4>  =r  Cri- 


-I  J?' 


, i c / -1,; 
V ft  "<f 


and  the  sides  are  (Cv£/  and 


The  slopes  of  the  diagonals  are  given  by 
1 


3 


U-iCQ 

v-Vb 


i/1 

f-i 

* J H s-«- 

ir=-ug 


*~L,  u--f~^J7L  l-tz 


7/ 


Lc  + ‘V*  t'o 


(28) 


In  this  case 
coefficients 
differential 
G ' dV2 . 


Ii!,  r> , and  G : are  identical  v:ith 
of  the  first,  fundamen tai0auadrat ic 
geometry,  dS2  - dX2  dXc  - E-dU2 


the 

form  of 
+ 2j? ' dUdV 


+ 
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and 


X 


-V~o 


k 


A = -U.0 
XT  - V"0 


■a/ 


C<-„'  l{-c~  */»  ^ 


The  major  diagonal  is  that  for  which  the  product  of  uQ  and 
vq  has  the  same  sign  as  F1.  Along  the  locus  of  F1  «=  0 the 
parallelograms  become  rectangles . 


If  u_  =»  v_,  the  lengths  squared  of  the  semi-diagonals 

O v# 

become 

*■-  [z ' ± V ■&')  u ’r 


and  their  slopes  are 


CL  ' 4-  CL. 
c//  / 


/ 

s. 


/ 

ci.,  — a 

2-/ 


( 

2.0- 


The  assumption  of  square  form  for  the  contour  P * 1,  In  the 
U#V  plane  13  not  Inconsistent  with  uQ  ^ vQ  because  we  may 
take  the  side  of  the  square  equal  to  the  longer  side  of  the 
rectangle-  All  the  points  then  lie  inside  the  square.  The 
only  loss  is  that  the  square  is  not  the  smallest  contour  for 
which  P * 1 . 


Graphical  representation  of  the  information 

The  assumption  of  bounded  errors,  while  crude,  may 
sometimes  be  near  enough  to  the  truth.  When  it  can  be  used 
together  with  the  condition  un  = v we  may  describe  the 
manner  in  which  the  transformation  propagates  errors  by  the 


-l  6- 


normalized  major  semi- diagonal 


/ 

z -+ 


Izf'i  ■*-  -y ' ) !/z" 


✓ 


which  is  the  quotient  of  the  longer  diagonal  of  the 
parallelogram  divided  by  the  side  of  the  square.  If  the 
contours  of  the  surface 


■f(. X,Y)  = £ <A  \Zf-'\  -h  & 


/ 


are  projected  on  the  X„Y  plane  they  may  be  regarded  as  loci 
of  points  for  which  the  circle  circumscribed  around  the 
parallelogram  is  of  constant  diameter.  In  other  words  they 
are  loci  of  constant  maximum  possible  error.  The  dis- 
advantage of  this  simple  representation  Is  that  a portion  of 
the  information  has  been  thrown  away.  We  have  discarded  the 
Information  about  the  two  directions  in  which,  for  uQ  - v , 
the  largest  measured  errors  may  yield  their  largest  and  smallesi 
displacements  in  the  X,Y  plane. 


A more  complete  way  of  presenting  the  facts  would  be  to 
plot  at  intervals  along  each  contour  of  the  normalized  major 
semi-diagonai , /<£f  -+  ) J?  *f')  + y j , a vector  showing  the 
direction  and  magnitude  of  the  normalized  minor  semi-diagonal 
and  a vector  showing  the  direction  of  the  major  semi-diagonal 

In  the  case  of  the  normal  la\«  of  errors  the  elements  of 
the  ellipses  for,  say  F - , can  be  represented  by  plotting 

small  L shaped  pairs  of  vectors  indicating  by  their  length 
and  directions  the  major  and  minor  axes  cf  the  ellipse:}  as 


they  vary  over  the  X,Y  plane. 


Measured  coord inates  not  of  the  same  physlca 1 dimensions 

So  far  In  the  discussion  we  have  supposed  that  U and  V 
are  of  the  same  kind.  If  they  are  not  of  the  3ame  kind  as, 
for  example,  in  ordinary  polar  coordinates,  there  are  no 
difficulties  in  applying  the  theory  for  normal  error 
distributions.  However,  since  the  moduli  of  precision  are 
then  of  different  dimensions,  it  is  futile  tc  discuss  a 
particular  case  In  which  k ~ 1c * . 

In  the  case  of  bounded  errors  the  preceding  treatment 
j.3  still  valid,  the  only  restriction  being  that  to  compare 
the  magnitudes  of  uQ  and  v^  is  without  meaning. 

Example 

For  the  sake  of  brevity  we  illustrate  the  theory  with 
an  example  of  such  simplicity  that  all  the  conclusions  are 
self-evident.  Consider  the  case  of  ordinary  polar  coordinates. 
Having  no  further  need  to  represent  the  differential  of  an 
error,  and  not  wishing  to  use  unfamiliar  symbols  for  the 
polar  coordinates,  we  write  the  transformation  as 


and  the  relations  between  the  errors,  written  as  d;i IToront.'.uls . 


as 


£ sz_  co-j  tP  ' ds  4-  A-*-+<  ■£  ' dt- 


*4--  - ~±^±.&k  a 


‘y 


The  equation  of  a normal  probability  contour  Is 


= ^-Jut-4-  idi-E?- 


— ! ' 4C~c~J~-S- 


"T-5 


si 


j Ax?- 4-  Z ' 2^l'(c - ~-l)  A:<id 


1~  Cl.  , 
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Hence 
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The  semi-axes  of  the  ellipse  are  therefore 

sJtu 

**  tf/ 

Note  that  the  former  is  independent  of  position;  and  the 
latter  is  proportional  to  r and  Independent  of  . 

The  orientation  of  the  ellipse  in  the  x,y  plane  is  given  by 
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In  other  words,  one  or  the  other  of  the  axes  of  the  ellipse 
is  always  parallel  to  the  radius  The  ellipses  become  circles 
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for  any  point  on  the  circle 


* = brJ 


Examining  the  sign  of  F,  we  find  that  inside  this  circle  F 
is  negative  in  the  first  and  third  quadrants  and  positive  in 
the  second  and  fourth.  Outside  the  circle,  F is  positive  :'n 
the  first  and  third,  and  negative  in  the  second  and  fourth 
quadrants . 

We  conclude  that  the  length  of  the  semi-axis  which  ;i.s 
in  the  radial  direction  is  uniform  and  equal  to  V<.  A-  • 

The  other  semi-axis  is  equal  to  /2L.  A"  ■ 

When  r - k’/k  the  ellipses  are  circle.-. 

We  next  replace  the  normal  law  of  errors  tilth  the 
assumption  of  bounded  errors  and  calculate  the  lengths  *:  i the 
diagonals  of  the  parallelograms  in  the  X,ir  plane  and 
inclinations.  The  lengths  of  the  semi-axes  are  given  by 


Since  there  Is  no  terra  in  drat,  the  pur-alicd.ogr-ar.jS  are 
everywhere  rectangles.  The  slopes  of  their  diagonals  a>.  $ 
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These  results,  in  such  an  extremely  simple  example, 
re  obvious. 


&crK. 

Donald  Poster 


